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ABSTRACT 



We study quantitatively the formation and radiation acceleration of electron- 
p ^, positron pair plasmoids produced by photon-photon collisions near Galactic black holes 

i (GBHs). The terminal ejecta velocity is found to be completely determined by the 

total disk luminosity, proton loading factor and disk size, with no dependence on the 
initial velocity. We discuss potential applications to the recently discovered Galactic 
superluminal sources GRS1915, GROJ1655 and possibly other GBHs. 
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Subject headings: acceleration of particles - accretion disks - gamma rays: observations 



INTRODUCTION 



o 

I The origin of relativistic jets from Active Galactic Nuclei (AGNs) remains one of the great 

challenges in astrophysics. Earlier Phinney (1982, 1987) concluded that radiation from the 
accretion disk is probably not responsible for accelerating the superluminal AGN jets with bulk 
^ ' Lorentz factors L ~ 10 (Porcas 1987). When F ^ 1, the aberrated photons in the jet's rest frame 



■ comes from the forward direction against the jet motion, thus limiting the maximum F achievable 

to only a few, much smaller than the observed values. 

Recently, observations of several Galactic black hole (GBH) candidates have revealed a causal 
relation between X-ray and radio flaring of the central object. In some cases, radio plasmoids 
were ejected relativistically away from the central sources following the x-ray flares ( [Mirabel et| 



al. 19921 , [Mirabel fc Rodriguez 199^ , [Hjellming fc Rupen 199^ . In two cases, GRS1915 ([Mirabel 



fe Rodriguez 1994[ ) and GROJ1655 ( [Hjellming fc Rupen 1995|) , the ejecta motions appear to be 



superluminal, resembling the superluminal AGN sources. 

Motivated by the correlation of the ejection with the x-ray flux, energetics of the ejecta and the 



relatively small F ~ 2.5, we proposed in a recent paper ( Liang &: Li 1995 ) a scenario for the origin 



of the episodic ejections based on the standard inverse-Compton accretion disk paradigm ( phapiro| 



et al. 1976). The idea is that episodic quenching of the (external and internal synchrotron) soft 



photon source to a level less than the internal bremsstrahlung flux superheats the innermost disk 
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to gamma ray temperatures ( [Liang &: Dermer 1988| , Wandel & Liang 1991, Pietrini & Krolik 



|1995[ ). However, these gamma rays cannot all escape, if their source compactness (oc gamma ray 
luminosity /source size, Svensson 1984, [Zdziarski 1984 ) is too high. Instead most gamma rays 
are converted into pairs via photon-photon pair creation which will be further accelerated axially 
by the radiation pressure of the disk x-ray flux and attain significant terminal Lorentz factors 
in a bipolar outflow if the disk luminosity is sufficiently high. When the ejecta encounters the 
interstellar medium and develops shocks and turbulence it will convert the bulk kinetic energy 
into high energy particles, emitting the observed radio synchrotron radiation. 

In this paper, we summarize some key results for pair plasma formation and their radiation 
acceleration by the disk photons for GBHs. 



2. COMPACTNESS AND FORMATION OF PAIR PLASMOIDS 



To study the formation of the pairs produced by 77 and 7X collisions outside the disk, we 
assume a 2-zone Keplerian disk model depicted in Fig.l. We define as the boundary, interior 
(exterior) to which local bremsstrahlung flux is higher (lower) than the soft photon flux. Gamma 
rays are emitted by a superheated ion torus within r-y ( [Rees et al. 1982 , Liang 1990| ), The spectrum 



from this region is assumed to be a Comptonized bremsstrahlung + annihilation spectrum with 
T-y ~ mcP (Fig. la, e.g. Zdziarski 1984, [Liang &: Dermer 198£ , Kusunose fc Takahara 198^ ). 



Exterior to r^, we assume that the disk x-ray flux has the standard inverse-Compton spectral 
shape (Fig. lb, Shapiro et al. 1976 , Sunyaev &: Titarchuk 1980| ). Using standard photon-photon 



pair production cross-section we compute the (angle-averaged) gamma ray pair-production depth 
function of photon energy due to collisions among gamma rays and between gamma rays 
and disk x-rays: 

r^^~r^ J d€'n{e)a^y{€,e)f (1) 
where a^^{e,e') is the angle-averaged pair production cross-section ( Gould &: Schreder 1967] ), and 



e = hu/bllkeV. n(e') = n^ie') + nx{e') is the normalized photon spectral density evaluated at r^, 
and / is a fudge factor of order unity to incorporate various complex geometry and angle-dependent 
effects (e.g. photon intensity anisotropy and variation with path length along the line of sight, 
angular dependence of cr^-y, etc.). We calibrate f using sample numerical simulations (details to 
be published elsewhere). For given geometry and spectral shapes (cf. Fig.l), r^^ basically depends 
only on the gamma ray source global compactness l^y = L^/r^/(3.7 x 10^^ erg/cm/s) ( pvensson 



1984[ , [Zdziarski 198^ ). In Fig. 2a we plot functions of photon energy for various l^ and 



~ mc^. In Fig. 2b we plot the fraction of gamma rays that are absorbed into pairs as a function 
of l^. We have assumed that all pairs produced escape without re- annihilation in situ due to their 
high kinetic energy (~ few hundred keV). This is a crude approximation but we estimate the error 
to be less than a factor of 2. The critical compactness at which ^ 80% of the gamma rays are 
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converted into escaping pairs is around 1^ ~ 100 (with ^ 2 uncertainty), consistent with similar 
results found by others (Zdziarski 1995, private communications). To achieve l-y > 100, Liang & 
Li (1995) show that the overall accretional luminosity Lq must be > 0.1 — 0.2 L^^- 



3. RADIATION ACCELERATION 



Consider a plasma of electrons, positrons and some entrained protons directly above the 
black hole and along the axis of the disk, interacting with the disk radiation via only Compton 
scattering (cf. Fig.l). Assuming that protons, electrons and positrons are always co-moving (a 



cold jet), 7p = 7e = r, and local charge neutrality is maintained (rig 
the equation of motion for the jet (at height Z) as: 



Up + ne+; 



we can obtain 



_ ^ 



( 



nip. 



m, 



eq 



) (^) Jdejcin iph{e,mi - - Pf^) - 1] + 



(2) 



GM/c^, and /x = cos 6, the angle between the photon momentum and the 



where z=k 

z-axis, and Iph{e,i^) is the spectral intensity seen by the jet (cf Fig.l, see P'Dell 1981] , phinneyj 



1987| for previous results). 



m. 



eq 



is the equivalent mass of the particles = rrif, + 



Hp /rip 
1, m, 



, the fraction of proton loading. For (^p 



0, m 



^-nir, with 



eq 



rUf 



for a pure pair jet and for 



eq 



m, 



+ nip for an electron-proton jet. In deriving equation (|^, we have: (1) kept 



Cp 

only the parallel (z) component due to the axial-symmetry of the system (though there will be 
a random orthogonal component); (2) used the angle- dependent differential Compton scattering 
cross section (in Thomson limit) and ignored the energy change (second order in scattered photon 
energy), but kept the momentum change (first order in scattered photon energy) for acceleration; 
(3) included the gravity pull on protons and leptons but neglected the radiation force on protons. 

The (isotropic) spectral intensity at the disk surface can be approximated as: 



Iphie) 



1.68 X 10^2 



1 1 



,-3 



e exp ( e) \LeddJ Mw Q 




1 



cm^ • s ■ sr 



(3) 



Here, 



where n = r/rg, Mio = M/IOM© and ^ = 18 - ^) + 2 (^Y^^ - (^Y 

V min max J \ max J \ min / 

rj^j„ and r^^j^ define the size of the disk from which the total luminosity is Lq, and we assume 



3/2' 



'^min ~ ^il'^g- have used the observed spectrum of GRS1915 ( [Harmon et al. 19941 ), but 
results are insensitive to detailed spectral shape since it will be integrated away when putting into 
equation (|2|) and only the total flux from each radius matters since we are working in the Thomson 
limit (as emphasized by the referee). The 7-ray component from r < (Fig. la) produces 
negligible additional effects on the axial acceleration due to the strong 7 — 7 self-absorption at 
high /-y. Substituting equation (^) into equation (^) the integral over becomes an integral over 
r,, (cf. Fig.l). 
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Equation @ can be evaluated numerically for the terminal Too- A key term in this equation is 
r^(l — — 1, which signifies the accelerating (or decelerating) force from a "ring" of radiation 
at radius r^,, depending on whether it is negative (or positive) for particles at the height z^, where 

= Zifj \J z1 + rl- One can find a disk radius r^ep, interior to which disk radiation accelerates 
the ejecta, and exterior to which radiation decelerates the ejecta. Whether an ejecta can get net 
acceleration depends on the sum of these two contributions. 

The dependence of Too on the initial conditions i^mitiai and zj^j„) and the geometry 
(^mm ^'^'^ "^max) ^au be obtained by solving equation We will focus on cases with Cp = first. 
The results are shown in Figure 3 and summarized here. 

(A) . For fixed disk geometry and total disk luminosity, regardless of their initial velocities 
^initiah leptous will be accelerated to a fixed terminal Too, even though particles with high T initial 
get decelerated first. The radiation field acts like a "thermostat" that regulates the particle's final 
velocity along the z-axis. In Fig. 3a, we show the ejecta's F as a function of height with different 
initial velocities, whereas fixing r^^^ = 12, r^^^ = 10^, z^^^ = 10 and Lq/L^.^^ = 0.2. 

(B) . Fig.3b, 3c, and 3d show the ejecta's F as a function of height with different rj^j„, 
different r^^^ and different initial z^i^, respectively. Again, Lo/Ledd = 0.2. It is evident that, Too 
depends somewhat sensitively on rj^^^ since it is mostly the inner disk radiation that accelerates 
the ejecta; on the other hand, the dependence on r^^^ is weak as long as it is » 100 due to the 
radiation intensity drops down drastically as r* increases. Note from equation (^) that the rate of 
acceleration is proportional to d^llph{e,il,), which is decreasing as 1/z'^. So efficient acceleration 
only occurs within a height of 100 — lOOOr^ of the black hole if the conditions are favorable. Fqo 
is not a sensitive function of zj^j^ as long as zj^j^ < 50. Ejecta starting at very high is not of 
interest here. 

Our main results are summarized in Figure 4, where we show the terminal Fqo as a function 
of Lo/Lf^dd for different proton loading (p = np/n^~ with rj^j„ = 12, r^^x = and z^i^ = 10. 
It is evident that for > 10^^, gravitational attraction on the protons hinders the acceleration. 
This has important implications for the composition of the ejecta. If the jets are composed of 
pure pairs as suggested in previous sections, then radiation acceleration can indeed account for 



the observed F in both GRS1915 and GROJ1655 ([Mirabel k Rodriguez 1994 [Hjellming fc Rupen 



19951) . The observed F ~ 2.55 of GRS1915 implies Lo/L^dd 0.2. This constrains the mass of 



black hole in GRS1915 to ~ I2M0 when using the observed luminosity ( Harmon et al. 1994 ) 



The jet could be suffocated if there is an additional isotropic radiation field liso and it reaches 
10% of 

Idisk- liso could be caused by the companion wind or circumstellar medium. 



4. DISCUSSIONS 



- 5 - 



The pairs are born with average kinetic energy of several hundred keV. However, we do not 
expect the "Compton rocket" effect ( P'Dell 1981 ) to be important to affect the ejecta's terminal 



velocity because the disk radiation field will regulate the particle's final velocity, no matter what 
the initial velocity is, as shown in Fig. 3a. 

We predict that the terminal Too should increase as L^/Ledd increases. This can be tested 
using the ejections from different episodes with the corresponding radiation luminosities. 

To summarize, in the present approach, we have adopted the thermal accretion disk model for 
the hard x-ray emissions (as suggested by the spectra of GRS1915, Harmon et al. 199^ ), and use 



that as the basis for the pair production and bulk acceleration. Hence the energetics of the jet is 
constrained by the pre-ejection accretional X-ray power output ([Liang fc Li 1995|) . We find that, 
for a jet composed mostly of leptons, disk radiation can accelerate it to a terminal Lorentz Too 
of a few, depending on the total disk luminosity and geometry. However, non-thermal processes 
may be needed to explain the power-law hard X-ray spectra of some x-ray novae (including 



GROJ1655, Harmon et al. 1994, Kroeger et al. 1995). Yet both classes of GBHs show relativistic 



outflows. Can GBH jets be also generated by say, purely electrodynamical processes unrelated to 
the conventional thermal hard x-ray disk models? We will explore alternative models in future 
publications. 

We thank an anonymous referee for constructive criticisms. This work is supported by NASA 
grant NAG 5-1547 and NRL contract No. N00014-94-P-2020. 
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Fig. 1. — Schematic depiction of the accretion disk along with emissions from the inner region 
(spectrum (a)) and the outer disk (spectrum (b)). At high compactness, spectrum (a) is not visible 
at infinity since most of it is converted into pairs within the dashed sphere. These pairs are then 
pushed out by the disk radiation with Lorentz factor T along the axis of the disk. 

Fig. 2. — (a) Pair-production optical depth function of photon energy. Curves a, b, and c 

are for compactness = 12, 100, and 392, respectively; (b) The fraction of gamma rays that are 
absorbed fabs as a function of l^. Note that fabs ~ 80% for = 100. 

Fig. 3. — Ejecta T as a function of height as solutions to equation (3) with Lo/Ledd = 0.2 and (^p = 
0.0 (pure pair jet). Parametric dependence of P on initial velocity ^initiah ^mim '^max a-i^d z^^^ are 
shown in plot (a), (b), (c) and (d), respectively. In each plot, other than the parameter that is 
varying, the rest are chosen from TiniUai = 1, r^,in = 12, r^^^ = 10^ and z*^-,^^ = 10. 

Fig. 4. — The ejecta's terminal T^o as a function of L^/L^dd for different proton loading C,p = np/rif.- . 
Curves a,b,c,d,e correspond to ("p = 0.0 (pure pair), 10^^, 10^^, 10^^, 1 (c — pjct). The insert 
corresponds to the vertical cuts for Lo/Ledd = 0.25,2.5 and 25, respectively. The observed P for 
GRS1915 is ~ 2.55. 
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